Threshold effects for two pathogens spreading on a network 
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Diseases spread through host populations over the networks of contacts between individuals, and a number 
of results about this process have been derived in recent years by exploiting connections between epidemic 
processes and bond percolation on networks. Here we investigate the case of two pathogens in a single popula- 
tion, which has been the subject of recent interest among epidemiologists. We demonstrate that two pathogens 
competing for the same hosts can both spread through a population only for intermediate values of the bond 
occupation probability that lie above the classic epidemic threshold and below a second higher value, which we 
call the coexistence threshold, corresponding to a distinct topological phase transition in networked systems. 



Social, technological, and biological networks of various 
kinds have been the subject of a large number of recent stud- 
ies published in the physics literature One of 
the principle practical applications of this body of work has 
been in modeling the spread of epidemic disease. Diseases 
spread over the networks of physical contacts between indi- 
viduals O, la 13 lal an d an understanding of the structure of 
these networks and the dynamics of disease upon them is cru- 
cial to the development of strategies for disease control. As it 
turns out, a large class of epidemicprocesses can be mapped 
onto bond percolation models allowing famil- 
iar techniques from statistical physics to be applied directly to 
their solution. 

An issue of some interest in current epidemiological re- 
search is the behavior of competing pathogens llll Il2l [l3i 
FA . Two diseases may compete for the same population of 
hosts because one disease kills hosts before the other can in- 
fect them. Or there may be cross-immunity between the dis- 
eases such that exposure to one disease leaves the host alive 
but immune to further infection by either disease. Compet- 
ing strains of influenza can show this type of behavior, for 
instance fl3i Q. The dynamics of competition between 
pathogens is in general complex, depending in particular on 
whether one pathogen gets a head start on the other in the 
population. In this paper we study the case in which two 
pathogens pass through the population at well separated in- 
tervals: one infects the population and causes an epidemic, 
leaving some fraction of the population immune or dead, and 
at some later time the second pathogen passes through the re- 
maining population. (Our arguments could also be applied 
to two successive outbreaks of the same disease.) The ques- 
tion we address is if and when the second disease is able to 
spread. If a sufficient number of hosts is removed from the 
population by the first disease then spread of the second be- 
comes impossible. As we will see there is a threshold value 
of the bond occupation probability or "transmissibility" for 
the first disease (a measure of contagiousness) at which this 
happens. This "coexistence threshold" coincides with a con- 
tinuous phase transition similar to the well known epidemic 
transition, but the two transitions are quite distinct: the co- 
existence threshold is an additional property of the network 
topology. 



Spread of both pathogens can occur only in the interme- 
diate regime between the epidemic and coexistence thresh- 
olds. Among other things, we determine by exact analytic 
calculation for a broad class of networks the position of the 
two thresholds. For the much-studied case of a "scale-free" 
network, we find that while the epidemic threshold for the 
first disease is always zero, the coexistence threshold is not. 
A corollary of this result is that while a single disease on a 
scale-free network cannot be eradicated solely by lowering 
the transmissibility, a similar intervention in the case of two 
competing diseases can eradicate one of the diseases, but not 
both. 

Consider then an epidemic taking place on a network of 
contacts between individuals. The network is represented by a 
graph in which vertices are individuals and (undirected) edges 
are contacts. The epidemic begins with a single individual 
and spreads along the contacts. Not every contact necessar- 
ily results in disease transmission however. We assume a 
generalized susceptible/infective/removed (SIR) dynamics for 
the disease of the kind described in |8] in which the disease 
spreads over edges with a probability T called the transmis- 
sibility of the disease. This dynamics can be mapped onto a 
bond percolation process on the same graph with bond occu- 
pation probability equal to the transmissibility 1^1 lil fioll . 
The connected clusters of vertices in the percolation process 
then correspond to the groups of individuals who would be 
infected by a disease outbreak starting with any individual 
within that cluster. Typically, for small values of T there are 
only small clusters and hence only small disease outbreaks. 
But above some critical transmissibility T c an extensive span- 
ning cluster or "giant component" appears, corresponding to 
an epidemic of the disease: once such a giant component is 
present, the pathogen reaching any of its members will infect 
them all and thereby reach an extensive fraction of the popula- 
tion. The value of T at which the giant component first forms 
is called the epidemic threshold and it corresponds precisely 
to the percolation threshold for percolation on the contact net- 
work. 

To be concrete, we examine in this paper the class of graphs 
that have specified degree distribution but which are other- 
wise random, in the limit of large graph size. (Recall that 
the "degree" of a vertex in a network is the number of edges 
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connected to that vertex.) Such graphs have been studied in 
the past by many authors 1151 lla Il7l 1 1 8t Il9ll and have be- 
come a standard arena for the exploration of epidemiological 
processes |7, 8,]. Epidemiological processes have also been 
studied on other types of networks, such as networks with de- 
gree correlations 1 20, 21], and it seems likely that the results 
presented in this paper could be generalized to such cases, al- 
though we do not do that here. 

Let pk be the fraction of vertices in our network that have 
degree k. We can also consider pk to be the probability that 
a randomly chosen vertex has degree k. The vertex at either 
end of a randomly chosen edge, on the other hand, has de- 
gree k with probability proportional not to pk but to kpk, the 
reason being that there are k times as many edges connected 
to a vertex of degree k than to a vertex of degree 1 , and hence 
the probability that our edge will be one of them is also multi- 
plied by k. We will primarily be interested in the distribution 
of the number of edges emerging from such a vertex other 
than the one we followed to get there. This excess degree is 
one less than the total degree of the vertex and therefore has a 
(correctly normalized) distribution 



qk = 



( k+ l)pk+i _ (k+ l)pk+i 



(1) 



where z = (k) = Y,k kpk is the mean degree of the vertices in 
the network. 

Our first pathogen can spread across the network if its basic 
reproductive number Rq is greater than unity, i.e., if for every 
person infected the mean number of additional people they 
infect is greater than 1 . 

When the disease arrives at a vertex, it has the chance to 
spread to any of the k other neighbors of that vertex, each of 
which chances is realized with probability T for an expected 
Tk additional vertices infected. Averaging over the distribu- 
tion qk of k, we find that the basic reproductive number is 
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Thus the disease spreads if and only if T is greater than the 
critical value 



T c = 



(k) 
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(3) 



To calculate the size of the epidemic when it does occur, 
it is convenient, following our previous approach I8t 1 1 71 12.211 . 
to define two probability generating functions for the distribu- 
tions pk and qk'. 



F o( x ) = £ PkX k , Fi (x) = qkx k = 
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where Fq denotes the first derivative of Fq with respect to its 
argument. In terms of these functions, for example, Eq. Q 
can be written 
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FIG. 1 : An epidemic of the first disease spreads through the network 
leaving some fraction of the vertices either immune to further infec- 
tion or dead, which we represent by simply removing them and their 
adjacent edges from the network (gray). The question we address is 
whether the remaining "residual network" (black) percolates and can 
therefore support the spread of the second disease. 



Now let u be the mean probability that a vertex is not in- 
fected by a specified neighboring vertex in the network during 
an epidemic outbreak of our disease. This quantity is equal to 
the probability that no transmission occurred between the two 
vertices, which is 1 — T, plus the probability that there was 
contact sufficient for transmission but that the neighboring 
vertex itself wasn't infected. The probability that the neigh- 
boring vertex wasn't infected is equal to the probability that it, 
in turn, failed to contract the infection from any of its k other 
neighbors, which is just u k with k distributed according to q^. 
Thus the mean probability of the neighbor being uninfected is 
YX=(>qk uk — F\ ( u )- Hence, u must satisfy the equation 



u = \-T + TFi{u). 



(6) 



Then the probability that a randomly chosen vertex is not in- 
fected is YX = QPk uk — Fq{u), and the fraction S of vertices that 
do get infected is one minus this: 



S=l-F (u) 



(7) 



Thus we can calculate the size of the epidemic by solving (|6jl 
for u and substituting the result into Q. S is also the proba- 
bility of an epidemic occurring if the disease starts with a ran- 
domly chosen individual — with probability 1 — S an outbreak 
fails to become an epidemic even when we are above the tran- 
sition. If S is regarded as an order parameter for the model, 
then the epidemic transition is a continuous phase transition 
in the mean-field universality class for percolation. 

Now consider the case in which our first disease causes an 
epidemic in the network, leaving a fraction S of vertices ei- 
ther dead or immune to infection by our second disease (or by 
a second wave of the first disease). To represent this mathe- 
matically, we remove these vertices from the network, leaving 
a smaller network of uninfected vertices which we call the 
residual graph — see Fig.^ Only if this residual graph has a 
giant component will it be possible for the second pathogen, 
provided it has a suitably high transmissibility, to spread. 

Clearly when T = for the first pathogen no individuals are 
infected and the entire graph remains for the second pathogen 
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to exploit. Conversely, when T = 1 an epidemic of the first 
pathogen will infect the entire giant component of the graph, 
and once this component is removed the second pathogen def- 
initely cannot spread (since, in the limit of large size, random 
networks have only one giant component). In between these 
two extremes, we can expect a transition, which we now inves- 
tigate. We begin by calculating the degree distribution of the 
residual graph. Once we have this distribution then, because 
the graph is uncorrelated, it is a straightforward exercise to 
determine whether it has a giant component or not. 

Consider a vertex with degree k. Let P(xmmf.,m\k) be the 
probability that it remains uninfected at the end of the first 
epidemic and has m edges that are attached to other uninfected 
vertices. In other words, P(uninf.,m|A:) is the probability that 
this vertex belongs to the residual graph and has degree m 
within that graph, given that it has degree k in the graph as a 
whole. 

This probability is equal to the probability that the ver- 
tex has k — m unoccupied edges that attach to infected ver- 
tices and m edges (occupied or not) that attach to uninfected 
vertices. The probability of an edge attaching to a unin- 
fected vertex is just F\{u) and the probability of being un- 
occupied and attaching to an infected vertex is (1 — T)(l — 
F\(u)) = u — F\(u), where we have used Eq. (j6ji. Then 
P(unmf. ,m\k) = Q [F x (u)] m [u - F\ {u)f- m . Multiplying by 
the probability p^ of having degree k and summing over k 
then gives the probability of being uninfected and having de- 
gree m within the graph of uninfected vertices: P(uninf.,m) = 

Zk=m Pk O [*i (")]'" [« ~ *l (")]*""■ Dividing by the prior 
probability P(uninf.) = 1 — S = Fq(u) of being uninfected, the 
probability distribution of the degrees of vertices within the 
residual graph is 



P(m|uninf.) = V p k ( ) [F x (u)] m [u - (u) 

F °( M ) k=m W 

The generating function for this distribution is 
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Given this generating function, we can determine whether 
the residual network has a giant component using the method 
ofRef. O. We define 



Fi(u) 



which is the generating function for the excess degree of a 
vertex reached by following an edge in the residual graph, 



precisely analogous to Eq. @. Then there is a giant com- 
ponent if and only if G\(l) > 1. Thus, the point G\(l) = 1 
constitutes an additional phase transition in the system, other 
than the standard epidemic transition, at which a sufficiently 
contagious second pathogen can cause an epidemic after the 
passage of the first through the network. We call this the co- 
existence transition and the point at which it occurs the co- 
existence threshold. Making use of Eq. (llOi . we find that the 
transmissibility T x at this point is the solution of the equation 



F 1 '(«) = l, 



(11) 



where u is a function of T via Eq. (|6}. 

For instance, in the case of a Poisson degree distribution 
for the original network p^ = e~ z z k jk\ (the standard Bernoulli 
random graph), we have Fq(x) = Fi (x) = e z (- v ~'), which means 
that the normal epidemic threshold falls at T c — l/z while the 
coexistence threshold falls at the point satisfying 



1 = z F l (u)=z(l-S). 



(12) 



If we can find S from Eq. (0, it is then a straightforward matter 
to find T x . 

The size C of the giant component in the residual graph, 
which sets an upper bound on the size of a possible second 
epidemic, is given by 

C=l-G (v), v = Gi(v), (13) 

as a fraction of the size of the residual graph Il7ll . To get 
the result as a fraction of the size of the original network, we 
then need to multiply by 1 — S. In Fig.|2]we show the sizes S 
and (1 — S)C of the epidemic and the giant component on the 
residual graph as a function of transmissibility for the Poisson 
case. As the transmissibility increases from zero, the size of 
the residual giant component is initially equal to the size of the 
giant component of the entire graph, which is very nearly 1 . 
As T passes the epidemic threshold for the first pathogen, 
however, the pathogen starts to spread and kills or renders 
immune to the second pathogen some fraction of the popu- 
lation, thereby reducing the size of the epidemic of the sec- 
ond pathogen. At some point — our coexistence threshold — so 
many are killed or made immune that too few are left to spread 
the second pathogen and C reaches zero. Thus the epidemic 
spread of both pathogens is possible only in the intermediate 
regime of transmissibility T c <T < T x indicated by the shaded 
area in the figure; if the transmissibility is either too low or 
too high, coexistence is impossible. In the inset of the figure 
we show how the two threshold values of the transmissibility, 
T c and T x , vary as a function of mean degree for the Poisson 
case. 

Of course, the mere existence of a giant component in 
the residual graph does not mean that the second pathogen 
will cause an epidemic. That depends on whether the trans- 
missibility of the second pathogen is high enough. Repeat- 
ing the analysis leading to Eq. (0, we find that the second 
pathogen can spread if its transmissibility is above the criti- 
cal value T' c = 1/Gj(l) or equivalently T'. = 1 /F'^u) — yet a 
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FIG. 2: The size of the epidemic of the first pathogen and the size 
of the residual giant component that it leaves behind, as a function 
of transmissibility on a graph with a Poisson degree distribution with 
mean degree z = 3. Inset: the position of the two thresholds as a 
function of mean degree for the Poisson case. The shaded areas in 
the plots denote the region in which both pathogens can spread. 



third threshold in our system (but one whose position is not 
solely a function of the network topology, since it depends 
also on the transmissibility of the first pathogen via Eq. 0). 
Noting that F[{x) is a polynomial with non-negative coeffi- 
cients and therefore monotonic increasing on the positive real 
line (within its radius of convergence), and that u < 1 since 
it is a probability, we see that F[(u) < F((l), and hence that 
T'c > T c - the minimum transmissibility necessary for the sec- 
ond pathogen to spread is never less than that necessary for the 
first. This accords with our intuition: as we have shown else- 
where [8], vertices with high degree are more likely to be in- 
fected than those with low degree, and therefore we would ex- 
pect the residual graph to have lower mean degree and hence 
higher epidemic threshold than the original network. 

Another example of interest is that of a network with a 
power-law degree distribution pk = k~ a /£(oc), for some con- 
stant a, where is the Riemann ^-function. Such networks 
are often called "scale-free." A variety of networks appear 
to be scale-free and they have attracted considerable atten- 
tion in the recent literature 0, • As is by now well under- 
stood 0,13, (uncorrelated) scale-free networks with a < 3 
have a vanishing epidemic threshold T c = because the sec- 
ond moment (k 2 ) of the degree distribution in Eq. Q diverges. 
Noting 1 17] that a power-law degree distribution gives a gen- 
erating function Fo(x) = Li a (.x)/^(a), where Li n (x) is the «th 
polylogarithm of x, and applying Eq. Jl li . we find by contrast 
that the coexistence threshold in such a network is in general 
nonzero. Furthermore, the critical transmissibility for the sec- 
ond pathogen T' c = \/F[(u) is also nonzero. (This nonzero 
threshold immediately implies that the residual network can- 
not itself be scale-free. The physical explanation of this result 
is that the first pathogen is more likely to infect higher-degree 
vertices and so selectively removes or immunizes the "hubs" 
in the network, destroying the power-law form. Removing 



hubs is well-known to be a good strategy for preventing the 
spread of disease 

The result T c = implies that a single disease spreading on 
a scale-free network of this kind can never be eradicated by 
an intervention whose sole effect is to reduce the transmissi- 
bility. Our findings indicate, however, that for the case of two 
competing pathogens on such a network, one of them can be 
eradicated by an intervention that lowers the transmissibility, 
but not both. 

To conclude, we have studied, using mappings to bond per- 
colation, the problem of two diseases spreading through the 
network of contacts between members of a host population. 
We find that, in the case where hosts can be infected with ei- 
ther one or other, but not both, of the diseases, the spread of 
both is possible only for intermediate values of the transmissi- 
bility of the first disease. There are two phase transitions that 
mark the boundaries of this intermediate regime. The first is 
the standard epidemic transition below which the first disease 
is not contagious enough to spread at all; the second is an ad- 
ditional topological phase transition in the network that cor- 
responds to the point at which the first disease removes from 
the population so large a fraction of the hosts that not enough 
remain to support the spread of the second disease. 

We have here studied only the simplest case of competing 
pathogens. A number of variants of the problem are of inter- 
est. For instance, in some cases the first pathogen may con- 
fer upon those it infects only partial cross-immunity to the 
second, so that the probability of infection with the second 
pathogen is reduced but not entirely eliminated. This pro- 
cess could be modeled using an extension of the formalism 
described here in which the residual graph is formed by re- 
moving a fixed fraction, randomly selected, of the vertices af- 
fected by the first epidemic. 
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